Spectroscopy, leptonic decays and the nature of heavy quarkonia 
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We examine the electronic width ratios of T resonances below the BB threshold by means of an effective (Cornell- 
type) QCD potential incorporating l/mb corrections obtained from a prior fit to the bottomonium spectrum. From 
our analysis we conclude that the T(25') and T(3S) states should belong to the strong-coupling (nonperturbative) 
regime while the T(IS') state should belong to the weak-coupling (perturbative) regime, in agreement with a 
previous study based on radiative decays. 



1. Introduction 

Heavy quarkonium has historically played a role 
of utmost importance in the rise of the Standard 
Model (SM), notably regarding the quark model 
of hadrons and the development of QCD as the 
presently accepted theory to describe the strong 
interaction among them. 

Furthermore, a large amount of data have been 
collected during the last decade at BEPC, B- 
factories, CESR, HERA, and Tevatron experi- 
ments, greatly improving the accuracy of the 
measured production cross sections, decay widths 
and branching fractions involving heavy quarko- 
nia (see [I] for a review). In the future, a Su- 
per Flavour Factory could provide further experi- 
mental results on heavy quarkonia to an unprece- 
dented accuracy |2l3j . 

On the other hand, such precise measurements 
are matched by the ever-growing soundness of the 
theoretical background, firmly based on an ef- 
fective field theory, namely the Non-Relativistic 
QCD (NRQCD). The maturity already reached in 
the field even makes feasible the search for new 
physics, e.g. in quarkonium decays, looking for ex- 
perimental deviations from the SM expectations. 
Let us mention the seek of light dark matter in in- 
visible quarkonium decays [4] (followed up by ex- 
perimental searches |5I6| ) and T radiative decays 
into dilcptons as a way of searching for a light 
non-standard Higgs boson [7I8I9I10I1I] . 
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With the advent of the quark model and QCD, 
hadronic properties have been traditionally un- 
derstood with the help of (more or less QCD- 
motivated) potential models, some of them having 
reached a fairly acceptable level in predicting or 
postdicting level spacings, transitions rates, etc. 
Nevertheless, potential models have several set- 
backs and limitations, mainly due to the fact that 
they do not come directly from first principles. At 
this stage, effective theories enter the game in or- 
der to describe rigorously the hadron dynamics. 

A prototype is the Heavy Quark Effective The- 
ory (HQET) which naturally describes hadrons 
with a single heavy quark |12| . These systems 
are characterized by two energy scales: the heavy 
quark mass, mg, and the characteristic scale of 
the strong interaction Aqcd. HQET is obtained 
by integrating out the scale rriq and expanding 
the QCD Lagrangian in powers of Aqcd. 

On the other hand, bound states made of two 
heavy quarks are characterized by more scales 
whose relevance and hierarchy are usually esti- 
mated by invoking the so-called velocity count- 
ing rules ^13j . Since the heavy quark relative ve- 
locity V is typically small (v^ ^ 0.3 for charmo- 
nium and ~ 0.1 for bottomonium) the different 
scales obey the useful relation i?„ ^ rnqv^ << 
p ~ rnqv << mg, where E„ is the heavy-quark 
bound-state energy with n the principal quan- 
tum number. NRQCD is obtained by integrating 
out the heavy quark mass uiq |14J. High-energy 
modes are not lost but encoded into short-distance 
coefficients and new local interaction terms in the 
effective Lagrangian. 
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The resulting framework allows the separation 
between the short-distance scale of the process un- 
der study from the longer distance scales associ- 
ated with quarkonium structure. Therefore opera- 
tors can be expanded as a double series (perturba- 
tive and nonperturbative), with and v being the 
expansion parameters controlling the accuracy of 
the truncated series, respectively. However, these 
expansion parameters are not completely indepen- 
dent in heavy quarkonium physics as the typical 
velocity of the heavy quark is determined by a 
balance between the kinetic energy raqv^ and the 
potential energy which should be dominated by a 
Coulombic-term ~ Ois/r. Setting r ~ 1/toqw, by 
invoking the virial theorem we are led to the well 
known relation 

asirnqv) ~ v 

Moreover, since as(/i) runs (decreasing) with 
the energy scale (for higher /i), one can approx- 
imately write 

asirrib) ~ 

where rUb denotes the bottom quark mass. The 
above relation should play an important role in 
assessing the velocity counting rules to be applied 
in quarkonium physics. In fact, relativistic cor- 
rections of order would be of the same order as 
perturbative corrections of order as{mi,). Let us 
also remark that the as perturbative expansions 
in NRQCD may be not fast convergent series and 
truncation even at NNL will likely imply sizeable 
effects [TS]. In this work, we will implicitly take 
the tts corrections into account through the veloc- 
ity counting rules. 

2. Potential model connection to pNRQCD 

A naive connection of NRQCD with potential 
models can be made by realizing that certain 
(long-distance, colour-singlet) matrix elements of 
NRQCD can be actually related in first approxi- 
mation to wave functions at the origin (WFO) or 
their derivatives. Nevertheless, this simple picture 
does not hold too far as there are NRQCD ma- 
trix elements without an equivalence in potential 
models, namely, colour-octet contributions since 
the heavy quark antiquark pair needs not to be 
in a colour-singlet state. Actually, the potential 
picture that emerges from NRQCD [TH] is quite 
different from the traditional one [18ll9j and su- 
perior. 



The observation that NRQCD still contains en- 
ergy scales irrelevant for the lower-lying states 
of quarkonium, led to further simplifications and 
the resulting theory was called potential NRQCD 
(pNRQCD) [20l21j . when only ultrasoft degrees 
of freedom remain dynamical. Such an effective 
field theory turns out to be, in fact, close to a 
Schrodinger-like description of the bound state 
|22I23| . Moreover, matrix elements in pNRQCD 
can be expressed as the product of WFOs and 
nonperturbative glue-dependent factors, yielding 
a formal similarity with potential models in many 
observable quantities (like decay widths into lep- 
tons or light hadrons [55] , or magnetic dipole tran- 
sition [H]), which do not happen in NRQCD. 

The relation between Aqcd and the scales rngv 
and mqv^ dictates the degrees of freedom of pN- 
RQCD: 

• The weak coupling regime, when rnqv^ > 
Aqcd > rnqv and the binding energy is 
mainly due to a Coulombic-like potential. 
Dynamics can be described using perturba- 
tive theory. 

• The strong coupling regime, when mqv ~ 
A-Qc D and the binding energy is mainly due 
to a confining (nonperturbative) potential. 

The assignment of each quarkonium state to any 
of these regimes is not such an easy task as the 
different scales are not directly measurable. The 
fact that the spectrum of excitations of the bot- 
tomonium family is not Coulombic suggests that 
the higher states are not in the weak coupling 
regime. However, there are claims in the litera- 
ture that the T(2S') and T(3S') can also be under- 
stood within the weak coupling regime [26 27 . In 
Ref.pH] T{1S,2S,3S) radiative decays were used 
to investigate their nature. Experimental results 
from CLEO were confronted with the theoreti- 
cal expectations, in particular as a ratio of decay 
widths for different energies of the photon. 

In this Letter we follow a similar strategy but 
focusing on electronic decays of the T(1S', 25", 35") 
resonances mainly basing our analysis on a phe- 
nomenological approach. Wc have employed the 
QQ-onia package presented in Ref. 29J to deter- 
mine the required WFOs, making use of a Cornell- 
type potential incorporating a -c'/r^ term as sug- 
gested by recent lattice studies. Table 1 shows 
the experimental values of the electronic partial 
widths for all three resonances. 
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Table 1 

Measured electronic widths (in kcV) for the T(IS'), T{2S) and T{3S) resonances (from |30|). 



TeeiTjlS)] r,e[T(25)] ree[T(35)] 

1.340 ±0.018 0.612 ±0.011 0.443 ± 0.008 



3. Velocity counting rules 

In this section we briefly review different ve- 
locity scaling rules appearing in both perturba- 
tive and nonperturbative regimes of heavy quarko- 
nia. All our arguments are order-of-magnitude es- 
timates which, moreover, are subject to uncertain- 
ties especially in the nonperturbative regime due 
to our ignorance of the scaling rules expected to 
be valid then. 

Notice that, no matter the regime, it always 
holds the following velocity counting rules 

{nS\ — \nS) - rribv^ ; {nS\V'^^^\nS) 
rub 

by the definition of v and from the virial theorem, 
respectively. The latter is an example where the 
naive guess of mtv is modified by the dynamics of 
the bound state to nibv^ . 

3.1. Perturbative regime 

The soft scale in heavy quarkonium is basically 
set by its size r which can be provided by the Bohr 
radius of the bound state. If r is small enough, i.e. 
r < I/Aqcd, the soft scale should be perturbative 
and the potentials can be entirely determined in 
perturbation theory. 

Then the following scaling rules should hold: 

{nS\V^^''>\nS) r^al/r^ r^mlv^ 
since as ^ v and r ~ l/(mbw), and 
{nS\V^^^\nS) - as/r^ - mgw'* 
and so on. 

Let us stress that perturbation theory cannot 
incorporate quark confinement, so it becomes cru- 
cial to determine the potential nonperturbatively 
in this regime. This can be the case for high- 
lying quarkonium resonances, as we are checking 
in this work. Indeed, states below the BB thresh- 
old and not too deep (namely T(2S') and T(3S')) 
are expected to be in the strong-coupling regime 
whereas the deeper T(IS') state is expected to be 
in the weak coupling regime. States above (or very 
close to) open bottom threshold are not expected 
to be in either regime [31| . 



3.2. Nonperturbative regime 

Admittedly, the power counting of NRQCD is 
not well known in the nonperturbative regime and, 
in fact, has been addressed by different authors 
in distinct ways [25132133] . We might assume a 
very conservative counting: nibv'^ with d standing 
for the operator dimensio Thus the following 
counting rules could hold: 

{nS\V^^^\nS) - mlv^ ; {nS\V^^'>\nS) - mfw^ 

On the other hand, making use of the l/rrib ex- 
pansion of the potential, one may write 

y(r) = y(o)(r) + ^-^ + — 4^ + ••■ (1) 

whereVW,fc = 0,l,2,.. . are the leading and sub- 
leading terms respectively. The potential T^'^^/m^ 
contains the leading-order spin-dependent poten- 
tials and the velocity-dependent potential. In this 
work we truncate the expansion up to y'-^^/m, ne- 
glecting 1/to^ terms and higher. In a nonpertur- 
bative regime this approximation should amount 
to a 0{v'^) accuracy at least when solving the 
Schrodinger equation, thereby justifying the use 
of a non-relativistic approach. 

Let us also point out that lattice calculations 
will be of help in determining the functional form 
of V^^^ /nib- We will come back to this important 
point in section 5.1. 

Finally, notice that the static potential V'^^''{r) 
is well parametrized by a Coulomb plus linear 
term (i.e. a Cornell- type functional form |35I36| 1. 

y(0)(r) = + crr + ^ 
r 

where a stands for the string tension governing 
the confining potential and ^ is a constant. This 
funnel shape will be used throughout this work 
as a reference, later on to be somewhat modified 
when defining the actual leading-order potential 
for heavy quarkonium. 

•^This is somewhat similar to HQET |34) where any opera- 
tor counts Uke Ai^^ 
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4. T leptonic decays 



5. Testing the nature of heavy quarkonium 



As is well known, leptonic partial widths are a 
probe of the compactness of the quarkonium sys- 
tem, and provide useful information complemen- 
tary to spectroscopy [T7]- In particular, the elec- 
tronic width for ^5';^ states {r[T{nS) e+e"]) 
probes the WFO according to potential models. 

Likewise, if the pNRQCD framework is applied, 
the leptonic width can still be written in terms 
of the radial WFO as mentioned in the Introduc- 
tion. The following expression obtained in Ref.pS] 
should hold up to order x {En/mi,, AQ(^jj/ml): 

Nc \Rnm^ 

2 ^ 

TT nit 



r[T(n5') ^ e+e"] = 



Im/ee("5i) 1 



En"^ 2e3 



2e: 



(2, EM) 



nib 

(0) 



3ml 



3ml 



mb 



(2) 



where Nc = 3 is a colour factor and En'' is the 
(leading-order) bound-state energy; e and B stand 
for universal (i.e. flavour and state independent) 
nonperturbative parameters, which can be ex- 
pressed in terms of gluonic field-strength corre- 
lators [25] . They can be determined either by ex- 
perimental data or by lattice simulation, but still 
their numerical values are quite uncertain. 

The matching coefficients Ivufeei^Si) and 
ImQeei^Si), corresponing to the Oii^Si) and 
Pi(^5i) operators of the NRQCD La grangian, are 
given at order by the expressions |14I37| 



Im/ee('5i) = i^Q'a' 



1 



IvCLQeeiC' Si) 



1 



OTT 



(3) 



(4) 



where \Q\ = 1/3 for the bottom quark. 

At lowest order, one recovers from (2) the well 
known formula |38j expressed in our notation as 

Nc Im/ee(^5i) \RI^\0)\' 



r[T(nS') 



e^e 



showing that the leptonic width of a quarkonium 
vector state is primarily sensitive to the square 
of its radial WFO, though perturbative and non- 
perturbative corrections are large indeed. Yet the 
unknown parameters in Eq. (2) do not allow its di- 
rect comparison with experiment. 



In the leptonic width ratio of two S-wave states, 
however, several terms cancel out, leading to 

T[TinS)^e+e-] |i?„(0)|2 



x[l 



(5) 



T[T{rS) ^e+e-] |i?r(0)P 

up to corrections of order 0{v'^), where q will be 
later determined in our analysis, providing an in- 
sight on the nature of heavy quarkonium and a 
hint at the velocity scaling rules to be applied. 
The correcting factor 5nr is given by 

fluigeei^Si 



2 esi^mb] 



Vlm/ee(35i) 9 
and using Eqs.(3-4) we get 

2^3(2mb)\ 
~^ 9 J 



^E^^^-Ei"^- 



miy 



E. 



(0) 



p(o) 



mb 



(6) 



where we have defined r^e up to 0{a^) corrections. 



1 



8Qs(2TOfc) 

3^ 



Now, taking into account that £"„ = M„ — 2mb 
with M„ being the meson mass [39] , we can safely 
use the following relation 

M"' ~ M"' ~ Er - E„ = Mr - Mr, 



Hence Eq.(6) can be rewritten as 
, 2e3(2mb) 





'Mr-Mn 




mb 



(7) 



On the other hand, tzi^) stands as the only 
nonperturbative gluonic parameter in (6), all oth- 
ers cancelling out in the ratio (5) at the desired 
accuracy. There is an experimental determination 
of this long-distance parameter in [1^ : 

£3(1 GeV) = l.StJ;? 

where the error bars are experimental only. Addi- 
tional theoretical uncertainties associated to sub- 
leading operators in the power counting and per- 
turbative expansion are not taken into account in 
the above uncertainty. 

Finally, we need to know 63 at the bottomonium 
scale. To this aim we can use the scale evolution 



e3(Ai') = e3(M) + 



2ACf ^^as{^l') 



(8) 



where Cp = 4/3, f3o = 11Ca/3 - 4n/TF/3 with 
Ca = 3, Tp — 1/2 and = 5 in our case. 

Setting 63(1 GeV) = 1.8, yields £3(2771,,) ~ 4.2 
which can be then used as an input in Eq.(7). 

^Let us note a factor Nq = 3 of difference between the 
definitions of the gluonic parameter e in 1401 and ea in 1251 . 
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Table 2 

Values of the predicted and experimental mass (in GeV), WFO squared (or derivative) in GeV'^+'^^, mean 
square radius (in fm) and typical quark velocity for the T(15', 25", 35*, 45), Xb(l-P:2P) and T{1D) states 
when the improved Cornell-type (Corn-mod) potential of Eq.(lO) is employed. For comparison, we present 
in the first column the masses obtained using a Cornell (Corn) potential [36] . 



Resonance 


Mass (Corn) 


Mass (Corn-mod) 


Exp. 




(r2)i/2 




T{1S) 


9.4603 


9.4603 


9.4603 


12.65 


0.23 


0.090 


Xb{lP) 


9.96 


9.8929 


9.9001 


1.409 


0.40 


0.071 


T{2S) 


10.05 


10.0236 


10.0233 


6.444 


0.51 


0.087 


T{1D) 


10.20 


10.1476 


10.1622 


0.562 


0.53 


0.078 


Xb{2P) 


10.31 


10.2729 


10.2600 


1.854 


0.63 


0.089 


T{3S) 


10.40 


10.3750 


10.3552 


5.404 


0.71 


0.103 


T(4S') 


10.67 


10.6477 


10.5794 


5.194 


0.88 


0.120 



5.1. Lattice estimates of yf^^/mt shape and 
size 

As commented previously, the V^^^ /mf, poten- 
tial has the form 1 /r^ relying on perturbation the- 
ory in the short-distance region. However, since 
the binding energy of the 6 — 6 system, typically 
of order mfoW^, can be similar or even smaller than 
Aqcd due to the non-relativistic nature of the sys- 
tem, it is essential to determine the potential non- 
perturbatively. Monte Carlo simulations of lattice 
QCD provide a powerful tool for a nonperturba- 
tive determination of the potential [41l42j . 

In a former analysis presented in Ref.[41j. 

was found to be comparable with 
the Coulombic term of the static potential (i.e. 
V^^^/mt, ~ l/r) when applied to bottomonium 
states up to r = 0.6 fm. Consequently, if F^^^ (r) 
is nonperturbative, the piece V^^'^{r)/m in the 
potential should not be considered as subleading 
with respect to y'^*'^(r). 

However, in a later study involving further long 
distance data up to r = 0.9 fm 42J, the same 
authors found that the l/r function was not sup- 
ported by the fit, while the functional form l/r^ 
with the linear term could fit the data well. It 
is interesting to note that l/^^^ /mt, turns out to 
have the same functional form as expected from 
perturbative theory. 

The y^^^ /m term cannot be neglected as com- 
pared to the static potential V^^^ and has to be in- 
corporated into the quarkonium potential for the 
sake of coherence. Thus the leading-order poten- 
tial should read 

1/(1) 

yLO = y(o) + i_ (9) 
nib 



yielding a potential of the form 
c c' 

Vcorn-rnodir) = :r + (7r + fj, (10) 

r 

We will refer to (10) as a Cornell-modified poten- 
tial, since the functional form has been improved 
by the additional —c'/r'^ piece; besides, the contri- 
bution from the V'^^^ jra term also alters the value 
of a [42]. 

In our approach the values of the parameters 
rribiC^d ,a and /i are obtained through a fitting 
procedure to the bottomonium spectrum (T(IS') 
and T(2S') states), not from lattice estimates. We 
obtain from the fit the following values for the 
parameters of the potential (10): 

a = 0.217 Gey^ c = 0.400, c' = 0.010 GeV"^ 

and TOfc = 4.7 GeVB 

The values of the predicted masses, WFOs and 
other properties of interest for different bottomo- 
nium states using this potential are shown in Table 
2. An excellent agreement with the experimental 
mass values of different resonances in the spec- 
trum can be observed. 

5.2. Discussion 

The WFOs corresponding to the Cornell- 
modified potential were obtained using our code 
based on a Numerov technique (See Ref.[29| for 
a thorough description of the QQ-onia package.). 
A comparison of the WFOs of the T( IS", 25, 35") 
states obtained using different potentials can be 
found in Table 3. 

significant coincidence is found between the mj, value 
used in the lattice calculation |42l and required in our fit. 
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Table 3 

|i?J^°(0)l^ (in GeV'^) values of T(ns) states (n = 1, 2, 3) for several potentials: BuchmuUer-Tye (mb — 4.88 
GeV) [11]; Cornell [18] {mt = 5.18 GeV); Cornell-modified (including a -c'/r^ piece, mi, = 4.7 GeV). 



Resonance 


|i?:;°(o)lB-T 




iRlf i^)\corn-mod 


T(IS') 


6.477 


14.07 


12.65 


T(25) 


3.234 


5.669 


6.444 


T(3S') 


2.474 


4.271 


5.404 



Let us remark that the inclusion of the 
—c'/r^ term in the Cornell potential has a non- 
straightforward effect on the new resulting WFOs, 
for it implies a modification of both the Coulom- 
bic term and the 6-quark mass obtained from the 
fit. Let us also note that the c' value is determined 
to a large extent by the T{2S) resonance in our 
fitting method (see [29] for more details). 

On the other hand, as already pointed out 
in section 5.1, the perturbative calculation of 
y(^) /rrif, yields the same functional r-dependence 
as suggested by lattice studies. Therefore, one 
may look upon V^^^-'/rrif, as an interpolating term 
between the perturbative and non-perturbative 
regimes, rendering the fit meaningful even using 
both T(1S', 25") states. As a check of the fit- 
ting procedure (see Table 2), the predicted me- 
son masses for the bottomonium family obtained 
from the Cornell-modified potential improve with 
respect to the Cornell potential [35], when com- 
pared with the experimental spectrum [30 . 

Moreover, the values of the WFOs shown in Ta- 
ble 3 look self-consistent: the (absolute) relative 
variations for the Cornell-modified versus the Cor- 
nell potential are ~ 10%, 14% and 26% for the 
T(15', 25', 35') resonances, respectively. Indeed, 
one could naively expect a 0{v'^) effect in the per- 
turbative regime as a consequence of incorporating 
the new term into the potential, but increasingly 
larger variations for higher (thus dominantly non- 
perturbative) states. 

Finally, note that the V'-^^/ml (and higher) 
terms neglected in the expansion (1) of the QCD 
potential should likely provide corrections of rela- 
tive order 0{v'^) to the WFOs obtained by solving 
the Schrodinger equation with the leading-order 
potential (10). This counting is supported by the 
spin-dependent splitting experimentally found in 
the bottomonium spectrum. In fact, if the static 
potential were exact, then the potential model 
would reproduce QCD up to corrections of rela- 
tive order v'^ HU . 



6. Numerical results 

Now we proceed to check the validity of formula 
(5) by recasting it onto the following double (ex- 
perimental to theoretical) ratio: 



r[T(n5) ^ ee]/r[T(r5) ^ ee] 

(|i?LO(0)|V|i?LO(0)|2)x[l + <S„, 



= 1 + A„, (11) 



where we have introduced the dimcnsionless quan- 
tity Anr {n,r = 1,2,3, n ^ r) parametrizing the 
deviation from unity for different combinations of 
aU three T(15), T(25), and T(35) states. 

The following experimental and theoretical in- 
puts have been employed in our analysis: 

1) The experimental input for the electronic 
widths can be readily obtained from Table 1, al- 
lowing a determination of the ratios with relative 
error ~ 1% 



r[T(n5) 



r[T(r5) 



(n, r = 1, 2, 3 n ^ r) 



2) The \R];°{0)\^ values for different potentials 
can be found in Table 3. We have assumed in (11) 
that 



|fl«(o)p 

|i?.(0)P 



where q is expected to be 2 on account of the 
arguments given in the previous section. 

3) 5nr were computed according to Eq.(7) (valid 
up to order (u^,^^) corrections). Experimental 
meson masses and m;, values for each potential 
were used (see Table 3); we set e3(2mb) = 4.2 
derived from Eq.(8) using 63(1 GeV) = 1.8 [40]. 

Thus, if the two quarkonium states n and r were 
in the same (strong) regime, one should expect 



100-w'J(%) 



(12) 
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Table 4 

A„r (in %) for different potentials from Table 3 using €3{2mb) = 4.2. Let us observe that A23 always 
remains smaller than A12 and A13, as expected if both states T(25', 35*) were in the same (strong-coupling) 
regime. 



Potential 


B-T 


Cornell 


Cornell- modified 


Ai3 


20% 


36% 


12% 


A12 


15% 


30% 


14% 


A23 


9% 


10% 


1% 



Conversely, one should expect A„r > 100 • v'^{%), 
if anyone of the states is in the strong-coupling 
regime and the other in the weak-coupling regime. 

The values for A„r obtained in our analysis for 
the BuchmuUer-Tye, Cornell and Cornell- modified 
potentials can be found in Table 4, representing 
our main result. 

In particular, it turns out that when the Cornell 
potential is employed we find that A 13 is greater 
than A12, which in turn is larger than A23, i.e. 

Ai3 ~ 35% > A12 ~ 30% > A23 ^ 10% 

in agreement with a counting rule providing q = 2 
for the latter case. 

In sum, our results for the Cornell potential are 
thus consistent with the expected level of accu- 
racy (up to order v^) of Eq.(ll), provided that 
both T(2S') and T(3S') states belong to the strong- 
coupling regime, while the T(15') state does not. 

Furthermore, once the — c'/r^ piece is included 
in the Cornell-modified potential we get 

Ai3^ A12 ~ 10% > A23 ^ 1% 

where A23 is now found to be remarkably small. 

Of course, there is an uncertainty coming from 
the 63 value in the computation of 5nr according 
to Eq.(7). Setting £3(2™?,) equal to zero, we find 
that A23 becomes appreciably worse. Thereby a 
non-null value of e3(2mf,) is clearly favoured in 
our analysis. Demanding A23 = in Eq.(ll) we 
get e3(2m6) = 3.2 from Eq.(7) and, consequently, 
€3(1 GeV) = 0.8 using the running equation (8). 

7. Conclusions 

In this Letter we have presented a phenomeno- 
logical study of the electronic width ratios of T 
resonances (below open bottom production), find- 
ing evidence favouring both T(25) and T(3S') 
states in the strong coupling regime, at the same 
time disfavouring the T(IS') in it, in accordance 



with the conclusions from the analysis of Ref. [28] 
based on radiative decays of T resonances. 

Moreover, the agreement between the experi- 
mental and predicted T(2S')/T(35') ratios is even 
better than naively expected from a conservative 
velocity counting once the Cornell potential be- 
comes improved by a —c'/r^ piece, motivated by 
recent lattice studies |41|42) . Therefore we can 
conclude that our results (both from lepton widths 
and spectroscopy) favour the inclusion of such a 
nonperturbative term into the Cornell potential. 

Let us also point out that a value of the gluonic 
nonperturbative parameter: £3(1 GeV) ~ 1.8^|!)'7 
(as found in [40j ) is compatible within errors with 
our analysis on leptonic decays. Actually, one 
might turn the question round extracting e3(2TOf,) 
from Eq.(ll) using the experimental data from Ta- 
ble 1 and the Cornell-modified potential, yielding 
e3(2m&) = 3.2 (£3(1 GeV) = 0.8) with an esti- 
mated uncertainty of ^ 30% assuming a ~ 10% 
accuracy in Eq.(7). Additional theoretical uncer- 
tainties should increase the allowed range though 
values of £3(1 GeV) around unity are preferred as 
a general result of our analysis. 

Finally, we want to stress the relevance of fur- 
ther accurate experimental measurements of lep- 
tonic widths (among other properties) of heavy 
quarkonia to carry out precise tests of effective 
theories of QCD (likely useful to deal with non- 
perturbative effects showing up at the LHC) and 
even direct searches for new physics |9,10 . A 
future Super Flavour Factory would play an in- 
valuable role in this regard. 
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